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ABSTRACT 

Understanding the origin and evolution of dwarf early-type galaxies remains an 
important open issue in modern astrophysics. Internal kinematics of a galaxy contains 
signatures of violent phenomena which may have occurred, e.g. mergers or tidal in- 
teractions, while stellar population keeps a fossil record of the star formation history, 
therefore studying connection between them becomes crucial for understanding galaxy 
evolution. Here, in the first paper of the series, we present the data on spatially re- 
solved stellar populations and internal kinematics for a large sample of dwarf elliptical 
(dE) and lenticular (dSO) galaxies in the Virgo cluster. We obtained radial veloci- 
ties, velocity dispersions, stellar ages and mctallicities out to 1-2 half-light radii by 
re-analysing already published long-slit and integral-field spectroscopic datasets using 
the NBuRSTS full spectral fitting technique. Surprisingly, bright representatives of the 
dE/dSO class {Mb = —18.0 ... — 16.0 mag) look very similar to intermediate-mass and 
giant lenticulars and ellipticals: (1) their nuclear regions often harbour young metal- 
rich stellar populations always associated with the drops in the velocity dispersion 
profiles; (2) metallicity gradients in the main discs/spheroids vary significantly from 
nearly flat profiles to —0.9 dex r~^, i.e. somewhat 3 times steeper than for typical 
bulges; (3) kinematically decoupled cores were discovered in 4 galaxies, including two 
with very little, if any, large scale rotation. These results suggest similarities in the 
evolutionary paths of dwarf and giant early-type galaxies and call for reconsidering 
the role of major mergers in the dE/dSO evolution. 

Key words: galaxies: dwarf - galaxies: elliptical and lenticular, cD - galaxies: evo- 
lution - galaxies: stellar content - galaxies: kinematics and dynamics 



1 INTRODUCTION 

Dwarf elliptical (dE) and lenticular (dSO) galaxies, low- 
luminosity {Mb > —18.0 mag) stellar systems with 
early-type morphology, represent the numerically dominant 
galaxy population in nea rby galaxy clusters and groups 
iFerguson fc Bingg:eli|[l993 ) . In the cold dark matter models 
of hierarchical galaxy formation they are considered as the 
building blocks of pr esently observed stellar systems (e.g. 
IWhite &i Frenkl Il99ll ). Thus, understanding their forma- 
tion and evolution becomes one of the important questions 
of modern astrophysics. Several dE/dSO fo rmation scenarii 
were proposed (see review and discussion in lDe Riicke et al] 
I2OO5I ). but none of them is able to fully explain all observa- 
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tional properties of these galaxies simultaneously. A possi- 
ble diversity of evolut ionary scenarii has been suggested by 
Ivan Zee et all (|2004ah and later investigated bv lLisker et al.l 
( 20081 ) to explain propertie s of different dE subclasses, al- 
though lBoselh et al.l (|2008al ) insist on the ram pressure strip- 
ping to be the unique channel of dE formation. 

Dwarf and giant early-type galaxies form two distinct 
sequences in the absolute magnitude { Mb) vs effective sur- 
face b rightness {{fJ-)e) diagram (see e.g. [Ferguson fc Binggelil 
I1994I ) joining around Mb ~ —18.0 mag, which is often 
referred as a lumin osit y separating these tw o classes (see 
iBoseUi et al]|2008bl and lKormendy et al.|[200i for the recent 
discussions). Although, so me arguments are provided (e.g. 
iGraham fc Guzmanll2003l ) for the continuity of the dwarf- 
giant sequence on the Mb vs (/ie) plot, it is clear that the 
effective surface brightness of dE galaxies correlates with 
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their luminosity, or, in other words, less luminous galax- 
ies are fainter in terms of surface brightness, making their 
observations a challenging task. One has to keep in mind 
that dE galaxies usually lack of the interstellar medium 
(ISM), therefore their spectra normally contain no emission 
lines. This explains wh y the first kinematical study of dEs 
(|Bender fc NietdllQQOl ') appeared more than a decade after 
the publication of the absorption-line kinematic s of giant 
elliptical galaxies (e.g. iBertola fc Capacciolilll975l ). 

New debates add r essing the dE formation appeared 
recently: iBoseUi et all (|2008al lbh argue for the ram pres- 
sure stripping to be the only scenario of dE formation, 
hence making these galaxies evolutionary different from gi- 
ant early-type systems an d reproducing t heir l ocus on the 
Mb vs plot; while IJanz fc Liskeij (|2008l ) succeeded 

to explain the observed discontinuity of giant and dwarf 
ellipticals in the size-luminosity relation by comparing a 
large homogeneous imaging dataset to the predictions of 
semi-analytical models of elliptical galaxy formation, thus 
suggesting their common origin. Comparing the structural 
properties of dwarf galaxies to their internal kinematics and 
stellar populations therefore becomes a question of major 
importance for understanding dE/dSO formation and evolu- 
tion. 

The significant progress in the astronomical instru- 
mentation and new detectors having low readout noise 
became the cornerstones to the accomplishments of sev- 
eral projects assessing spatially-resolved kinematics of dif- 
fuse elliptical galaxies in different envir onmental condi- 
tions. Kinematical profi l es presented in De Riicke et al.l 
ll200ll,l2003l):lGeha et al] (l2002l.l2003l):|Pedraz et all (|2002D : 
ISimien fc Prugniell (|2o'ol ): Ivan Zee et al.1 (|2004bh revealed 
a diversity of the degree of rotational support and pres- 
ence of k incmatically-decouplcd component s in some ob- 
jects llpe^ Riicke ct al...2004. : Gcha ct al...2005l : IPrugniel et~al] 

Detailed studies of the mor- 
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I2OO5I : iThomas et al] 
phological properties of dEs revealed embedded structures 
such as stellar discs, sometimes harbo uring low-contrast 
spiral arms, or bar s in many ga l axies (|jerien et al.l I2OO0I : 
iBarazza et al] |2002| : iLisker et al.1 l2006h . These discoveries 
strengthen the evolutionary connection between dwarf el- 
lipticals and discy late-type dwarfs. 

At the same time, studies of the stellar population prop- 
erties of dwarf elliptical galaxies still remain a challenge both 
for observations and for the data analysis, because absorp- 
tion line strengths, a tradit ional spectroscopic i ndicator of 
stellar populations (see e.g. IWorthev et al.|[l994l for the de- 
scription of the Lick system) require high signal-to-noise ra- 
tios in order to reach reasonable quality of age and metallic- 
ity estimates. Some attempts were made to deri ve dE stellar 
popu lations from the narrow-band photometry (|Rakos et al.l 
120011 ). However, the calibrations used by the authors were 
based on Galactic globular clusters, therefore an assumption 
about very old ages was done, resulting in low metallicity es- 
timates (—1.4 to —0.6 dex). The integrated measurements 
of dE stellar populations from the Lick indic es were made 
later (|Geha et al.ll2003l : Ivan Zee et al.l 12004^ 1 and revealed 
populations having intermediate ages of 3-5 Gyr, and con- 
siderably higher metallicities around —0.5 dex. These results 
were confirmed and extended to a larger sample of Virgo 
cluster dEs recently (|Michielsen et al.llioosi '). 

A new family of techniques for recovering stellar 



population prope rties from integrated ligh t spectra ap- 
peared recently (ICid Fernandes et al.l l2005l : lOcvirk et al.l 
I2OO6I : IChilingarian et al.1 12OO7X " These methods are based 
on the direct fitting of the the stellar population models 
against the observed spectra in the pixel space. Since every 
pixel of an absorption-line galactic spectrum bears informa- 
tion about the stellar content of a galaxy, these techniques 
are potentially much more sensitive than the line-strength 
indices dealing with the individual spectral features. An- 
other advantage of the pixel space fitting is an easy way of 
getting rid of residual cosmic ray hits and regions of the 
spectra contaminated by emission lines by excluding these 
regions from the fitting procedure. 

The full spectral fitting using the NBurst s code 
(see preliminary descri ption in IChihngarianI l2006l : 
IChilingarian et al.1 bOOTbl ) allowing to extract internal 
kinematics simultaneously with the parametrized star 
formation history, was successfully used in the studies of 
integrated kinematical and st ellar population prop erties 
of dE gala xies in the Fornax (iMichielsen et all |2007| ) and 
Abell 496 (|Chilingarian et al. I l2008d ) clusters. The method 
was shown to be stable and working well even for the data 
having low signal-to-noise ratios. 

Up-to now, a published sample of dwarf and low- 
luminosity early-type galaxies with the spatially-resolved 
stellar population information (apart from Local Group 
dwarfs studied by means of colour-magnitude diagrams) 
contains only seven objects, all o f them observed with the 
Multi-Pupil Field Spectrograph (|Afanasiev et al.l lioOll ) at 
the 6-m Bolshoi Teleskop Al'tazimutal'mj (BTA) of the Spe- 
cial Astrophysical Observatory of the Russian Academy of 
Sciences (SAO RAS ). The sample includes f our dE/dSOs 
in the Virgo cluster ([Chilingarian et al]|2007^[d[) and three 
low-l uminosity E/SOs in groups: NGC 770 (IPrug niel ct al.l 
l2005h . NGC 126 and NGC 130 jChihngarian et a l. 2008^ 
In all cases the two-dimensional maps of kinematics and stel- 
lar population properties were obtained using the NBuRSTS 
technique. Chemically and evolutionary decoupled struc- 
tures which a re quite common in giant elliptical and l entic- 
ular galaxies (|SirchenkQ.2006l : iKuntschner et al.ll2006h were 
detected in the central regions of six low-luminosity galaxies. 

Presently, the growing amount of high-quality optical 
spectroscopic data on dwarf galaxies becomes available in 
the archives of major observatories. Therefore, we decided 
to re-analyse the existing and available spectroscopic data 
for dE/dSO galaxies using the full spectral fitting technique. 
The project aims at studying the connection between inter- 
nal kinematics and stellar populations of dE galaxies and 
understanding the mechanisms of their formation and evo- 
lution, primarily in the cluster environment. 

In the first paper of the series we provide the kinemat- 
ical and stellar population profiles of Virgo cluster dE/dSO 
galaxies obtained from the new analysis of several published 
datasets carrie d out using the NBursts full spectral fit- 
ting technique ([Chilingarian et "al]l2007bl l3). The thorough 
discussion of the obtained results will be given in the next 
papers. The paper is organized as follows: in Section 2 we 
describe the spectroscopic data, data reduction procedures, 
and define a sample of galaxies. Section 3 provides details 
about the techniques used to analyse the data. Section 4 
presents the results of this analysis, the conclusions are given 
in Section 5. 
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2 THE DATA AND THE SAMPLE 

In the course of our study we have used intermediate- 
resolution ground-based spectroscopic data coming from 
various sources. The criteria we applied to select the datasets 
were: (1) only published data with available advanced cali- 
brations such as stellar templates, twilight spectra, or well- 
defined spectral resolution information and well documented 
spectrograph; (2) coverage of the blue/green spectral region 
(4800A< A <5400A) containing prominent absorption fea- 
tures; (3) intermediate spectral resolution not lower than 
R = 1300 to be able to constraint kinematics simultaneously 
with stellar populations and tackle the degeneracy between 
metallicity and velocity dispersion; (4) we gave preference 
to the data having intermediate or high signal-to-noise ratio 
and spatially resolved datasets (i.e. longslit or IFU spec- 
troscopy). Hereafter we present the data collections used in 
order of their importance for our study and briefly discuss 
essential data reduction steps which had to be considered 
prior to the analysis. 

2.1 Palomar dE Project: Long-Slit Spectroscopy 

Observations of 16 Virgo Cluster dwarf and low-luminosity 
early-type galaxies were conducted with the Double Spec- 
trograph at the 5-m Palomar telescope during two observ- 
ing runs in March 2001 and April 2002. The reduced, sky 
subtracted flux calibrated datasets were kindly provided by 
L. Van Zee (PI of the project) in fall 2004. The two wave- 
length ranges were observed with the double beam spectro- 
graph: 4800-5700A with the resolving power R « 2200, and 
8250-8900A having higher spectral resolution {R « 5400). 
Integration times ranged from 1200 sec for VCC 1075 to 
12000 sec for VCC 1308. The 2 arcsec wide slit spanned 
2 arcmin and was place d along major axes o f the observed 
galaxies (see Table 2 in Ivan Zee et al.|[2004bl ). A number of 
radial velocity standards (F7-K2 giants) was observed. 

We refer to al l the details r e garding observations and 
data reduction to Ivan Zee et al.l (|2004bl ). where thorough 
kinematical analysis of these data is presented. In this study 
we make a direct use only of the blue beam spectra cen- 
tered on the Mgfc tr iplet; the red beam kinematics from 
Ivan Zee et al.l (|2004bl ) is used for comparison with the kine- 
matical profiles we obtain. 

We used the spectra of HD 136202, HD 182572, and 
HD 187691 obtained with the same instrument in March 
and September 2001 and provided to us together with the 
spectra of galaxies, to compute the variations of the spectral 
line spread function (LSF) along the wavelength and assess 
the precision of the data reduction. The analysis of the 10 
spectra of HD 182572 obtained for different positions of the 
star on the slit allowed us to assess the spatial variations of 
the spectral resolution. 

We fit the stellar spectra against the high-resolution 
spectra of the corresponding stars available in the 
ELODIE.3.1 (|Prugniel et al.ll2007l ) library in five wavelength 
segments covering the spectral range of the blue setup of 
the Double Spectrograph and overlapping by 20 per cent. 
The fitt ing was done using the penah zed pixel fitting pro- 
cedure (|Cappellari fc EmsellemI l2004h returning v, cr, /is, 
and ^4 coefficients of the Ga uss-Hermite parametrization 
llvan der Marel fc Franxl [l993h . The heliocentric radial ve- 
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Figure 1. Wavelength-dependent variations of the Gauss- 
Hermite parametrization of the Double Spectrograph LSF at the 
Palomar 5-m Hale Telescope obtained from the fitting of the 10 
HD 182572 spectra. The adopted curves are displayed by solid 
bold lines. 



locities available from the ELODIE.3.1 library have been 
subtracted from the fitting results. Then the heliocentric 
corrections were computed for every individual spectrum us- 
ing the IRAF NOAO.RV.RVCORRECT task and applied to the 
measurements. 

In case of perfect and unbiased wavelength calibration 
this procedure would result in zero radial velocities, while 
non-zero values would measure the systemic errors of the 
dispersion relation. Variations of a, hs, and /14 along the 
wavelength are caused by the distortions introduced by the 
optical system of the spectrograph both, in the collima- 
tor and in the camera. The mean value of the /14 coeffi- 
cient anti-correlates with the slit width: wider slit result- 
ing in the H-shaped line profile c orresponds to negative ^4 , 
as in e.g. VLT FLAMES/GiraflFe (|Chilingarian et aLll2008^ ). 
where the physical diameters of optical fibers are larger than 
the diffraction limit of the collimator. 

In Fig [1] we present the variations of the Gauss-Hermite 
parametrization coefficients along the wavelength obtained 
from the fitting of 10 individual spectra of HD 182572. The 
adopted behaviour of the coefficients which we use to trans- 
form the PEGASE.HR models into the resolution of the 
Double Spectrograph are shown by solid bold lines. 

The measurements for individual spectra are scattered 
between —10 and 30 km s^^, which is probably caused by 
the fiexures in the instrument resulting in the overall shift of 
the wavelength solution zero-point. Therefore, we can expect 
systematic errors of the radial velocity measurements of an 
order of 20 km s^^. However, since our study is not aimed 
at the accurate measurements of absolute radial velocities of 
galaxies, this effect is not critical for us and affects none of 
our scientific conclusions. We did not notice any correlation 
between the measured systemic radial velocity offset and the 
position of the star on the slit. Good agreement of the radial 
velocity profiles for some objects with different datasets (see 
below) also suggests no important systematic errors of the 
wavelength solution along the slit. 

We also see an overall trend along the wavelength range, 
which is very well reproduced by the measurements made on 
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all individual spectra: the systemic values first smoothly de- 
crease by ~7 km s~^ between 4800 and 5400A, and then 
increase by ~8 km s~^ to the red end of the wavelength 
range. This behaviour may be a result of an insufficient or- 
der of the polynomial used to approximate the wavelength 
solution during the data reduction procedure. 

The spectral resolution (cjjjg^) as a function of wave- 
length (Fig[T]) rapidly rises from ~65 to ~55 km s^^ between 
4800 and 5050A then slowly increasing to ~48 km s~^ at 
5650A. This behaviour keeps very similar for all the analysed 
spectra, suggesting no important variations of the spectral 
resolution along the slit. 

Within the precision of the measurements (~0.03), the 
hs coefficient remains stable along the wavelength range 
with the mean value of -1-0.05, suggesting slight overall asym- 
metry of the DS blue arm LSF. The /i4 coefHcient also 
demonstrates very stable behaviour for different wavelength 
with the mean value of 0.00 with a smooth gradient from 
—0.03 to -f 0.03, which is the same order as the uncertain- 
ties of the measurements. 



Hence, the overall variations of the Double Spectro- 
graph LSF in the blue beam exhibit very reproducible be- 
haviour along the wavelength showing neither significant 
changes along the slit, nor time-dependent drifts. The wave- 
length dependent variations of the spectral resolution and 
radial velocity systematics are taken into account while fit- 
ting the spectra of galaxies by convolving the PEGASE.HR 
stellar population models using the wa velength-dependent 
kernel, as described in Section 4.1 of IChilingarian et ahl 
l|2007d ). Briefly, several (usually, five) sets of models are cre- 
ated by convolving the original grid with LSF for different 
(five) positions along the wavelength, then the linear inter- 
polation between these sets is done at every wavelength. 
This allows to obtain trustable velocity dispersion measure- 
ments down to ~l/3 of the instrumental resolution (i.e., 
~20 km s~^) at intermediate signal-to-noise ratios (~20) 
which is very important fo r reliable estimations of th e stel- 
lar population parameters (jChilingarian et "al]|2008al ). 

The initial sample of Ivan Zee et all (|2004bh included 
two galaxies having poor signal-to-noise ratios: VCC 1743 
and VCC 1857. For these objects we provide only integrated 
measurements of kinematics and stellar populations along 
the slit. For the other 14 objects we have performed the 
adaptive binning along the slit in order to reach the min- 
imal signaLto-noise ratio of 20 or 30 (10 for VCC 1075) 
per individual bin per pixel at 5200A. This was done by 
co-adding individual consequent pixels along the slit from 
the centre towards outer regions of the galaxy until the tar- 
get signal-to-noise ratio had been reached. The first bin was 
always assigned to the photometric galaxy centre on the 
slit and grown in both directions. For the kinematical and 
stellar population analysis the data have been rebinned log- 
arithmically in wavelength with the step corresponding to 
40 km s-^ 



2.2 HyperLeda FITS Archive: Long-Slit 
Spectroscopy 

The HyperLeda projeclQ (|Paturel et al.ll2003l ') known mostly 
by its largest existing database of homogenised structural, 
photometric, and kinematical parameters of galaxies, also 
contains a rich collection of raw and processed spectral and 
imaging datasets known as the HyperLeda FITS Archive 
(HFA). This resource includes all the data on early- type 
galaxies presented in a series of papers by Simien & Prug- 
niel originating from the CARELEC long-slit spectrograph 
at the 1.93 m telescope of I'Observatoire de Haute Provence. 
In our study we use the subset of this collection devoted to 
dwarf and low-luminosity early- type galaxies in th e Virg o 
cluster and presented in detail in lSimien fc Prugniell (|2002h . 
The data were obtained during 5 observing runs in April and 
June, 1999, February and March, 2000, and January, 2001. 
The CARELEC setup #2 provided intermediate spectral 
resolution of i? « 5000 at the wavelength range 4700 < A < 
5600A for the 5 arcmin long slit. The slit widths as speci- 
fied in the metadata (FITS headers) were: 2.2, 2.1, 2.3, 1.8, 
and 1.7 arcsec for the five observing runs mentioned above. 
The integration times ranged from 60 minutes in a single 
exposure for NGC 4476 to 210 minutes in four exposures for 
I C 3461. The complete obse rving log is available in Table 2 
of ISimien" fc Prugniell (|2002l ). 

We used the original raw CCD frames available in the 
archive and reduced them using our software package origi- 
nally developed for the IFU data reduction and adapted to 
the long-slit spectroscopy. The main steps of the data reduc- 
tion include: (1) subtraction of bias; (2) cleaning c osmic ray 
hits using Laplacian filtering (|van Dokkum|[200ll ) (3) pre- 
liminary flat-field correction using dome flat; (4) automatic 
identification of the arc line spectrum (He, Ne, Ar) and 

rd 

building the wavelength solution using the 3 order polyno- 
mial; (5) rebinning the spectra logarithmically in wavelength 
with the step of 25 km s^^; (6) residual flat fleld correction 
using twilight frames. In average, 30 arc lines were used to 
build the wavelength solution, resulting in the accuracy of 
about 0.04 A. 

For all the observing runs, but April 1999, we used the 
twilight spectra obtained in the same observing mode of the 
instrument in order to assess variations of the CARELEC 
LSF along the wavelength range and along the slit. We fit 
them against a high-resolution spectrum of the Sun available 
in the ELODIE.3.1 spectral library exactly the same way as 
the stellar spectra from the Palomar dE project. The LSF 
was shown not to exhibit any changes along the slit. At 
the same time there are modest variations along the wave- 
length range: the radial velocity systemic offset and fjjjg^ 
(in parenthesis) change smoothly from +2 (32 to 50 for dif- 
ferent observing runs) in the blue to —2 (25 to 40) km s~^ 
in the red end of the wavelength range; stays within 0.01 
from zero; /i4 decreases from -f0.04 to 0.00. 

There are important changes of the instrumental res- 
olution o-jjjg-)- between the observing runs because of dif- 
ferent slit widths used. No twilight spectra were obtained 
in the April 1999 observing run, however, two stars in- 
cluded in the ELODIE.3.1 library were observed: HD 84937 
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and HD 137759. The former is a relatively early type star 
(F5 III), therefore its spectrum does not contain sufficient 
amount of deep spectral lines to measure accurately the 
spectral resolution variations. However, the latter one, a K- 
type giant, has an appropriate spectrum to determine the 
spectrograph's LSF. 

We did not use the data for four galaxies, IC 3120, 
IC 3457, IC 3653, and NGC 4482. IC 3120 is a late-type 
dwarf and is out of the scope of this study; IC 3653 was 
observed under very poor atmosphere conditions (6 arcsec 
seeing); and the two remaining objects have too low signal 
(perhaps due to poor transparency for NGC 4482) to de- 
termine stellar population parameters. We notice that the 
data for IC 3457 also contain the spectrum of PGC 41516 
(VCC 1399) on the slit, but it is also too faint even for in- 
tegrated measurements of stellar populations. 

Finally, we used the data for 15 galaxies, 3 of those 
also having the data available from the Palomar dE project. 
Spectra of some objects (e.g. IC 783A) have poor signal-to- 
noise ratios, therefore we provide only integrated measure- 
ments of the stellar population parameters along the slit. 
For the remaining galaxies we performed the adaptive bin- 
ning along the slit to reach the minimal signal-to-noise ratio 
of 10 or 15 per bin per pixel at 5200A using the technique 
described in the previous subsection. 

2.3 IFU Spectroscopy of Four dE Galaxies 

Four bright dE/dSO galaxies, IC 783, IC 3468, IC 3509, 
and IC 3653 were observed using the MPFS IFU spec- 
trograph at the Russian 6-m telescope BTA SAO RAS 
during 3 observing runs in March 2004, May 2004, and 
May 2005. The details about obs ervations, da t a redu ction 
and analysis were presen ted in IChiUngarianI (|2006h and 
IChilingarian et al. I (l2007dldh. All the data are available in 
the ASPID archive (jChilingarian et al.|[2007al l. 

Here we include the original measurements presented 
the corresponding papers for two galaxies, VCC 490 and 
VCC 1871. We re-analysed the data for another two ob- 
jects, VCC 1422 and VCC 1545, improving the sky subtrac- 
tion and increasing the power of the multiplicative poly- 
nomial continuum to 15 (see next section). The data were 
spatial ly rebinned using the Voro noi adaptive binning tech- 
nique (|Cappellari fc CopinI [20o3) to reach minimal target 
signal-to-noise ratio per bin. 



2.4 SDSS DR6: Aperture Spectroscopy 

The Virgo cluster is partially contained in the foot- 
print of the Sloan Digital Sk y Survey Data Release 6 
(|Adelman-McCarthv et al.ll200l l. Therefore, relatively high 
signal-to-noise ratio flux calibrated intermediate-resolution 
[R ~ 1800) spectra covering a large wavelength range 
(3800 < A < 9200A) are available for many dE/dSO clus- 
ter members. We use SDSS spectra for the objects from 
the samples presented above in order to compare the stellar 
population properties in the inner regions of the galaxies. 
Relatively low spectral resolution of the data did not allow 
us to measure velocity dispersions below 40 km s^^. 

The spectra were obtained using the 2.5 m Apache Point 
Observatory telescope with the multi-object double spec- 
trograph containing 640 3-arcsec wide circular fibres in the 



3 deg field of view (|Gunn et al.ll2006h . The spectroscopic 
target selection of SDSS is based on the limiting magnitude 
inside the fibre aperture, therefore the sample is biased to- 
wards nucleated dwarf galaxies. 

Since neither twilight spectra were available in the 
SDSS archive, nor bright stars included in the ELODIE.3.1 
library were targeted by the SDSS due to high fluxes, it 
is not possible to directly measure the LSF variations of 
the SDSS spectrograph along the wavelength. However, the 
changes of the spectral resolution in terms of o-jj^g-j- are com- 
puted for each fibre by the SDSS data processing pipeline 
and provided together with the data (in the 5^^^ extension 
of FITS files containing extracted spectra), fjjjg^ decreases 
from 68 km s"^ at 3800A to 52 km s"^ at 5800A. Then it 
sharply jumps to 66 km s^^ smoothly decreasing down to 
51 km s~^ at 9200A. This happens because SDSS spectra in- 
clude two parts (blue and red) obtained in the two different 
arms of the instrument. Here we fit the SDSS spectra in the 
wavelength range between 3900 < A < 6700A (rest-frame). 
In addition to a wide coverage of the spectral domain re- 
sulting in better quality of the stellar population parameter 
estimates by the full spectral fitting, the centre of the se- 
lected wavelength interval, 5300A, well corresponds to the 
discussed long-slit spectroscopic data. 

23 of 31 galaxies included in our final sample have SDSS 
DR6 spectra available. 

Our final sample is presented in Table [T] We provide 
the common galaxy designations, their n umbers in the Virgo 
Cluster Catalogue (|Binggeli et al.|[l985^ and availability of 
data in every of the discussed datasets. The sample spans 
a range of luminosities —18.27 < Mb < —15.10 mag, cen- 
tral velocity dispersions 23 < a < 79 km s^^. It is not a 
statistically representative sample of the Virgo cluster dE 
galaxies, because of the well-known correlation of the lumi- 
nosity and effective surface brightness of dE/dSOs, making 
fainter galaxies (in terms of the luminosity) more difficult 
to observe. For several objects among the least luminous 
members of our sample we provide only integrated mea- 
surements of the stellar populations properties. The faintest 
object with the spatially resolved information available is 
VCC 1308 {Mb = -15.41 mag). 



3 DATA PROCESSING AND ANALYSIS 



We h ave fit the high-resolution pegase.hr (|Le Borgne et al.l 
l2004h simple stellar population (SSP) models against the 
observatio nal data using the NBurs ts full spectral fitting 
technique (|Chilingarian et "aP '2007c''l?). We computed the 
SSP models for the lSalpeteii (| 19551 ) stellar initial mass func- 
tion us ing the recent versio n of the ELODIE stellar library 
(v 3.1. |Prugniel et al.ll2007^ providing a spectral resolution 
of i? = 10000 in the wavelength range 3900 < A < 6800 A. 

The spectral fitting method is based on the non- 
linear least square fitting using the constrained Levenberg- 
Marquardt optimisation of the following quantity: 



2 _sr^{F,- Pip{T,{SFH) (g) C{v, (7, ha, hi) + P2,))' 



where Ti{SFH)= ^ kiT,{t^,Z„) (1) 
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Table 1 . Final sample of low-luminosity dwarf early- type galax- 
ies in the Virgo cluster. Columns (1) and (2) give the object name 
and VCC number, columns (3) - (5) indicate availability of the 
data, column (6) provides the source of photometric data. Galax- 
ies observed with the MPFS IFU spectrograph are shown in italic. 



Name 


VCC 


Pal. 


OHP 


SDSS 


Ph.ref* 


IC 3081 


178 


yes 


no 


yes 


3 


IC 781 


389 


no 


yes 


yes 


4 


UGC 7399A 


437 


yes 


no 


no 


1 


IC 783 


490 


no 


yesl 


yes 


4 


UGC 7436 


543 


yes 


yes 


yes 


1 


IC 783A 


545 


no 


yes^ 


yes 


4 


NGC 4328 


634 


no 


yes 


yes 


4 


IC 3328 


856 


no 


yes 


yes 


1 


IC 3344 


917 


yes 


no 


yes 


2,3 


IC 3363 


965 


yes 


no 


yes 


3 


IC 3369 


990 


yes 


no 


yes 


3 


NGC 4431 


1010 


no 


yes 


no 


4 


NGC 4436 


1036 


yes 


yes 


yes 


2,3 


IC 3383 


1075 


yes 


no 


yes 


1 


IC 3381 


1087 


no 


VPS 


VPS 


1 


IC 3393 


1122 


yes 


yes 


yes 


3 


IC 3413 


1183 


no 


yes 


no 


4 


NGC 4476 


1250 


no 


yes 


yes 


1 


NGC 4482 


1261 


yes 


no^ 


yes 


1 


IC 3437 


1308 


yes 


no 


no 


3 


IC 3461 


1407 


no 


yes^ 


yes 


1 


IC 34 68 


1422 


no 


no 


yes 


1 


IC 3486 


1491 


no 


yes 


yes 


4 


PGC 41726 


1514 


yes 


no 


no 


3 


IC 3509 


1545 


no 


no 


yes 


1 


IC 3602 


1743 


yesl 


no 


yes 


1 


IC 3647 


1857 


yes^ 


no 


no 


1 


IC 3653 


1871 


no 


no'^ 


yes 


1 


IC 3735 


2019 


yes 


no 


no 


1 


IC 3773 


2048 


no 


yes 


no 


1 


IC 3779 


2050 


yes 


no 


yes 


1 



^ poor signal-to-noise ratio, only integrated measurements are 
possible 

^ very poor signal-to-noise ratio, object excluded 

^ poor seeing, object excluded 

* Photometric refe rences: 1 — [ F errares e et al ] ll2006h. 2 - 
IStiavelli et"al] 1I2OOII) and/or lGeha et al.l ||2003|) , 3 - Ivan Zee et al.l 
(|2004al l. 4 - SPSS DR6 



where £. is the line-of-sight veloc ity distribution in the 
Gauss-Hermite parametrization (|van der Marel fc Franxl 
Fi and AFi axe observed flux and its uncertainty; 
Ti{SFH) is the flux from a synthetic spectrum, represented 
by a linear combination of Ntursts SSPs and convolved ac- 
cording to the line-spread function of the spectrograph; Pip 
and P2q are multiplicative and additive Legendre polyno- 
mials of orders p and q for correcting the continuum; t„ 
and Z„ are age and metallicity of the n*'' SSP. The rela- 
tive contributions of the SSPs and the additive continuum 
are fit linearly at every evaluation of the function, whereas 
all the remaining parameters (including P29 for the additive 
continuum) are derived from the the non-linear optimisa- 
tion. The needs for multiplicative continuum and possible 
side-effects have been presented and analys ed in detail in 
A ppendix A2.3 of IChilingarian et al.1 (l2007d'l. Append ix Bl 



If emission lines present in the spectrum, the wavelength 
regions contaminated by them are excluded from the fitting. 
In the wavelength range of SDSS we then exclude Balmer 
lines (7A-wide intervals), [OIII] (A = 4363, 4959, 5007 A), 
[Nil] ( A = 5199 A), [NI] (A = 6548, 6583 A), and [SII] (A = 
6717, 6731 A). Besides, we always exclude from the fitting 
the spectral regions contaminated by the strong air-glow 
lines: [OI] (A = 5577, 6300, 6364 A), Nal (A = 5890, 5896A), 
and mercury lines originating from the hght pollution, the 
strongest one being Hgl (A — 5461 A). 

The stellar population properties were derived in the 
parameter space with the r otated age and me t allicity axes 
as described in Section 3 of IChilingarian et al.l (|2008j ): 



A ppendix J\z.6 otiuniimgarian et al.l ll/UUml. Append i 
of IChilingarian et all (|2008cl ). and lKoleva et all (|2008l ) 



r? = (3Z-f21ogiof)/V13; 

6» = (-2Z + 31ogl(,^)/^/T3 (2) 

In this coordinate system the rj axis is parallel to the di- 
rection of the well-known age-metallicity degeneracy for old 
stellar populati ons in the age-metallicity space expressed by 
IWorthevI (|l994l ) as "if the percentage change Aage/AZ ~ 
3/2 for two populations, they will appear almost identical 
in most indices". 

All the results presented in this paper correspond to 
the fitting of a single SSP. We did not include the additive 
polynomial in the fitting. The adopted values of the multi- 
plicative polynomial continuum were 13, 13, and 15 for the 
Palomar dE, HFA, and SDSS datasets respectively. In one 
of the next papers of the series we will present the fitting 
of the nuclear regions of the galaxies with a more complex 
model containing tw o stellar popu l ations . 

It was shown bv lKoleva et al.l (|2008l ) that the full spec- 
tral fitting using the NBuRSTS code produced consistent 
stellar population parameters with those derived from the 
measurements of the Lick indices, being several times more 
precise. Here we also stress that (1) age determinations from 
the full spectral fitting are insen sitive to the masking of 
Balmer lines (see A ppendix A2 in Chilingarian et al.ll2007cl 
and Appendix B in IChilingarian et al.ll2008cl ): (2) velocity 
dispersions can be precisely determined at down to 1/3-1/2 
of the spectral resolution (i.e. 15-20 km s^^ for the HFA 
data and 18-25 km s~^ for the Palomar dE data) and these 
measurements re main unbiased (see e.g. [ Kpleva et al.ll2007l 
and Section 2 in IChilingarian et al.1 l2008ah : (3) more op- 
timal usage of the information contained in the spectrum 
by the full spectral fitting compared to the techniques us- 
ing line-strength indices results in considerably higher pre- 
cision of the derived stellar population parameters or, recip- 
rocally, lower values of the required signal-to-n oise ratio per 
A to get the comparable uncertainties (see e.g. iKoleva et al] 
[2003); (4) luminosity- weighted values of ages and metallic- 
ities are insensitiv e to the q/Fc rati os of the popul ations 
being fit (fchilinga rian et al.|[200 8c: Koleva et al.ll2008l ) even 
when using the PEGASE.HR SSP models not represen- 
tative of metal-rich a-enhanced populations since they are 
based on the empirical stellar library, which contains stars 
only from the Solar neighbourhood, where a/Fe abundance 
ratios are correlated with their metallicities. 

We have conducted a numerical experiment to quan- 
tify the age- and metallicity-sensitive information in the 
absorption- line spectra in order to illustrate the points (1) 
and (3) mentioned above. Since the fiux at every wavelength 
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contributes to the total value of and this contribution 
depends on many parameters, a possible way to estimate 
the sensitivity of the fitting procedure S(A,po, • • ■ ,Pn) to a 
given parameter pi at a given wavelength A is to compute 
the corresponding partial derivative of the template grid ac- 
counting for the multiplicative continuum variations. In our 
case the sensitivity to a stellar population parameter r) (and 
similarly, 6) at a point of the parameter space {to, Zo,ao) 
will be expressed as; 



s{\-n)\{to,zo.cTo) 
d 



t=5 Gyr, [Fe/H]=-0.3 dex, (7=60 km/s 



x'{PiA\n:e)\(,to,zo){n\,t,z)®c{ao))}, (3) 



with the quantities defined in Eq. [T] [5] assuming SFH con- 
taining a single SSP. This means, that the squared fitting 
residuals of a spectrum (r;o + A77, 6^0) against the model 
at every pixel would reflect its contribution to the 
overall when varying rj, i.e. its sensitivity to, or, in other 
words, the measure of information about a given stellar pop- 
ulation parameter. 

Practically, the computation is straightforward and 
has been conducted as follows. We have chosen six PE- 
GASE.HR "reference SSPs" representing stellar popula- 
tions with ages of 1.5, 5, and 12 Gyr and metallicities of -1.0 
and -0.3 dex. Then, for every of the six "reference SSPs", 
six models have been constructed by convolving them with 
three Gaussians corresponding to the velocity dispersions of 
60, 120, and 180 km s^^ and adjusting rj and 6 parame- 
ters corresponding to their ages and metallicities (see Eq. [5)) 
by 0.05 dex. For example, for the (5 Gyr, —0.3 dex) "ref- 
erence SSP" the pairs were constructed using (5.33 Gyr, 
—0.258 dex) and (5.50 Gyr, —0.328 dex) respectively. Later, 
these models were fit a gainst their "reference SSPs" using 
the PPXF procedure bv lCappellari fc EmsellemI (|2004l 'l with 
the 13th order multiplicative polynomial continuum. The fit- 
ting residuals obtained by this procedure correspond to the 
quantity defined in Eq. |3l 

We have co-added the information on r)-, and 6- 
sensitivity in 20 A bins in the wavelength range between 
4700 and 5600 A and normalised it by the total value of non- 
reduced 'x^ thus obtaining the relative importance of every 
20 A-wide bin to the determination of the stellar popula- 
tion parameters. Four examples with different stellar popu- 
lations and velocity dispersions are provided in Fig [2] One 
should not directly compare blue and red curves, because 
they represent the normalised quantities: the absolute val- 
ues of 5(A, rf) as defined by Eq. |3] are several times higher 
than ^(A, 6') explaining the elongated shapes of the 1-a un- 
certainty ellipses in the age-metallicity space. 

As expected, the most valuable information is contained 
in prominent absorption-line features such as H/3 and Mg6 
having corresponding Lick indices defined. However, it is no- 
ticeable that rj and 0-sensitive information is present every- 
where in the spectrum although in different amounts. Go- 
ing to low metallicities or young ages, thus, reducing the 
absorption-line features apart from H/3, causes concentra- 
tion of information in the most prominent lines. Increasing 
the intrinsic velocity dispersion also does so by smearing out 
faint spectral features, therefore for high-resolution spectra 
of low velocity dispersion objects it is possible to take the full 












2logt + 3[Fe/H] 
3logt - 2[Fe/H] 















t=5 Gyr, [Fe/H]=-1 .0 dex, (7=60 km/s 




t=5 Gyr, [Fe/H]=-0.3 dex, cr=180 km/s 




Wavelength, A 

Figure 2. Spectral distribution of the stellar population sensitive 
information for four stellar populations in 20 A-wide bins. Cor- 
responding age, metallicity, and velocity dispersion are indicated 
in every panel. The studied SSPs are shown in black; blue and 
red histograms display the relative importance of spectral bins 
(in per cent) for the determination of -q and 8 stellar population 
properties respectively. Green bars denote the bands defining the 
following Lick indices: H/3, Mg6, Fe527o, and Fe5335 with yellow 
side-bars corresponding to the corresponding pseudo-continuum 
definition regions. 



advantage of techniques such as NBuRSTS over line-strength 
indices. 

To be stressed that in none of the presented examples 
the 20 A-wide bin around H/3 contains even 20 per cent of 
the ^-sensitive information. This explains why it is possible 
to precisely constrain the stellar population age using full 
spectral fitting te chniques even if H/3 (or oth er Balmer lines) 
are not available (|Chilingarian et al]|2008al H). 

To measure the metallicity gradients we fit the co-added 
spectra in the two regions along the radius adding the data 
on both sides of the galaxy centre, individually selected for 
every object. With our data it is not possible to define the re- 
gions automatically based only on the values of r^, because 
the effective surface brightness values as well as exposure 
times differ quite significantly from object to object, leading 
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to very different signal-to-noise ratios at a given distance 
from the centre in terms of r^. Therefore, we based the def- 
initions of the radial zones based on the observed radius of 
a nucleus r„uc defined as a half-size of a pronounced bump 
in the metallicity profile, if observed or 1.5 arcsec other- 
wise, resulting in a range of values from 1.5 to 4 arcsec (for 
VCC 478). The inner radius of the inner region was cho- 
sen to be ru = 1.5r„uc to prevent the contamination of the 
spectra by the light of a galactic nucleus, but not smaller 
than 3.5 arcsec. The outer radius of the inner region and 
the parameters of the outer region were defined as follows: 
rio = (2.5/1.5)rii; Toi = (3. 5/1. 5)rii; Too = (5.5/1.5)rii. 
At the end, inner and outer regions turned to be about 
1/4 to 1/2 and 1/2 to 1 in the units of effective radii. 
The values of the metallicity gradients were computed as 
A[Fe/H]/Ar = 2([Fe/H]„ - [Fe/H]i)/(ro, + Too - [n^ + no)), 
where [Fe/H]i and [Fe/H]o are the values of metallicity de- 
rived from the fitting of the co-added spectra in the inner 
and outer regions along the radius. Then, the values were 
converted into AFe/H per and per kpc. 

We computed the radial variations of the B- 
band stellar mass-to-light rat ios using the pegase.2 
iFioc fc Rocca-Volmerangell997l 'l code for the Salpeter IMF. 
It is important to mention that according the the predictions 
of PEGASE.2, in addition to the stellar content, another 30 
to 40 per cent of the mass is supposed to be contained in 
the form of ISM. The cluster dwarf elliptical galaxies are 
supposed to efficiently lose their ISM in quit e short time 
by m eans of gravitational harassment (s ee e.g. Moore et al.l 
1 19961 ) or ram-pressure stripping (see e.g. lAbadi et al.lll999[ r 
Therefore, to estimate the mass-to-light ratios of the stel- 
lar populations we do not take into account this gas mass. 
The two other remaining components of the total SSP mass, 
white dwarfs and neutron stars / black holes are taken into 
account. 

The parameters referred throughout the paper as "pe- 
riphery" or "main galactic body" were obtained from the 
co-added spectra along the slit on both sides of the photo- 
metric nuclei of galaxies at radii from ra to roo- 

For the objects having spectra with sufficient signal- 
to-noise ratios in the blue spectral range we also measured 
the Lick indices. We compared the measurements of Mg6 
and the combined iro n index (Fe) = 0.72Fe527o + 0.28Fe5335 
(|Thomas et al.ll2003h with the SSP models for different a- 
enhancement presented in the same study. These measure- 
ments were made in the nuclear and peripheral parts of the 
galaxies, as well as for all the SDSS DR6 spectra. To measure 
the Lick indices we have transformed the data into the Lick 
system by convolving the spectra with the quadratic differ- 
ence of the Lick spectral resolution and the effective spectral 
resolution of the data (instrumental resolution and the in- 
trinsic broadening d ue to velo city dispersion of the galaxies) , 
as was proposed by Kuntschncr (20041') and applied to the 
SDSS data, e.g. in IChilingar"ian et*^ (|2008d 'l. PI We derived 
the values of the a/Fe abundance ratios by inverting the 



^ Since we did not have sufficient number of observations of 
Lick standard stars, no systemic offsets between our data and 
the original Lick/IDS response function (see e.g. Appendix A of 



grids of models presented in iThomas et al] (|2003l ') using the 
age values obtained from the full spectral fitting procedure 
described above. 



4 RESULTS 

In Table [2] we present the kinematical and stellar popula- 
tion properties of the galaxies in our sample. We provide 
the global photometrical parameters, absolute magnitudes 
Mb corrected for the Galactic extinction and half-light radii 
?-£ from the literature. The conversion from the photomet- 
ric systems in the original bibliographic so urces to the B 
band was done using the transformations of lFukugita et ahl 
1 19951 ). We compiled the photometric measurements from 
the followin g studies (see Table [T ]): the ACS Virgo Clus- 
ter Survev |Ferrarese et al.l |2006| ). the Keck dE project 
' Geha et al.ll2q02l.l2003ll ^T WFPC2 surfa ce photometry 



Stiavelh et all |2001^ . P alomar dE project (Ivan Zee et al.l 



iMichielsen et all |2008|) could be calculated and applied to our 

measurements of the Lick indices 



l2004bl P). SDSS DR6 (|Adelman-McCarthv et al.l |2008| ). 

From the SDSS DR6 catalogue we have adopted the values 
of the p'-band 50-percentile Petrosian radius as r^- We have 
used the distance modulus of the Virgo cluster —31.09 mag, 
as in lFerrarese et al] (|2006l ). and respectively corrected the 
absolute magnitudes from other studies. 

The following quantities are provided for the nuclear 
regions and main bodies of the galaxies: velocity dispersion, 
luminosity-weighted age, metallicity, stellar mass-to- light ra- 
tios assuming the Salpeter IMF, [Mg/Fe] abundance ratio. 
We also provide metallicity gradients per kpc and r^, com- 
puted as described above. 

Since the stellar population properties estimated by the 
full spectral fitting are a subject to the age-met allicity de- 
generacy, we provide the uncertainties in the rotated coordi- 
nate system (see Eq. [5]). The two pairs of columns, A?; and 
A6 provide a possibility to recover the approximate dimen- 
sions of the inclined error ellipses in the age-metallicity plot 
for the nuclear and peripheral parts of the galaxies. 

4.1 Comparison of the Datasets 

For the long-slit and IFU datasets we have fit the co-added 
spectra of the central regions of the galaxies in the 3 arcsec- 
wide aperture in order to compare the derived stellar popula- 
tion properties directly with those obtained from the fitting 
of the SDSS spectra. In Fig [3] we compare these measure- 
ments. The agreement is very good showing no systematic 
offsets between the different datasets, which suggests that 
the data were reduced properly. Since there is no system- 
atic offset in the metallicity estimates, we conclude that the 
subtraction of additive background components, such as dif- 
fuse light in the spectrograph or the sky background, was 
precise. This is an important conclusion, that SDSS spectra 
allow us to study the stellar population properties of dE/dSO 
nuclei. Although the spectra are contaminated by the light 
of main galactic discs/spheroids, the results obtained from 
the full spectral fitting of SDSS spectra are fairly consistent 
with those obtained from the analysis of other ground-based 
data, giving possibility to study much broader samples of 
dwarf early-type galaxies using only SDSS data. 

For the three objects, VCC 543, VCC 1036, and 
VCC 1122 having long-slit data in both, Palomar dE project 
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Table 2. Kinematical and stellar population properties of dwarf early-type galaxies: circumnuclear regions. The columns are as fol- 
lows: object name, absolute magnitude (extinction corrected), effective radius, adopted nuclear radius, central velocity dispersion, age, 
metallicity, 2 columns representing uncertainties of stellar population parameters (see text), central [Mg/Fe] abundance ratio, B-band 
mass-to-light ratio, 2 columns providing ages and metallicities from the fitting of SDSS DR6 spectra. 

VCC Mb re r„ ao t [Fe/H] Ar? AS [Mg/Fe] {M/L)b isDSS [Fc/H]gDss 

mag arcsec arcsec km s ~^ Gyr dcx dcx dex dex {M/L)q Gyr dcx 



0178 


-15.43 


8.7 


2.0 


50±2 


10.6 


-0.73 


0389" 


-16.28 


10.9 


2.5 


23±3 


3.9 


-0.38 


0389° 


-16.28 


10.9 


2.5 


36±2 


4.9 


-0.30 


0437 


-16.85 


25.1 


4.0 


47±1 


4.5 


-0.36 


0490" 


-16.14 


17.4 


17.4 


33±5 


5.1 


-0.23 


0490'' 


-16.14 


17.4 


1.5 


34±8 


3.3 


-0.35 


0543 


-16.55 


17.1 


2.0 


49±1 


5.4 


-0.49 


0543" 


-16.55 


17.9 


2.0 


48±3 


5.7 


-0.41 


0545" 


-15.10 


9.5 


9.5 


31±7 


5.1 


-0.47 


0634" 


-16.31 


16.6 


2.5 


26±2 


4.0 


-0.30 


0856" 


-16.42 


14.7 


2.5 


31±1 


5.0 


-0.43 


0917 


-15.48 


7.8 


2.0 


41±1 


4.4 


-0.62 


0965 


-15.44 


13.8 


1.5 


34±4 


2.2 


-0.56 


0990 


-16.31 


9.7 


2.0 


53±1 


4.0 


-0.40 


1010" 


-16.79 


22.9 


2.5 


43±2 


7.3 


-0.33 


1036 


-17.03 


14.6 


1.5 


53±1 


3.5 


-0.15 


1036" 


-17.03 


14.6 


1.5 


41±2 


2.6 


-0.02 


1075 


-15.71 


14.6 


1.5 


53±7 


5.7 


-0.74 


1087" 


-16.80 


18.9 


2.5 


63±4 


5.3 


-0.27 


1122 


-16.03 


11.4 


2.5 


40±1 


2.4 


-0.37 


1122" 


-16.03 


11.4 


2.5 


33±1 


2.4 


-0.30 


1183" 


-16.66 


10.5 


2.5 


56±2 


4.7 


-0.28 


1250" 


-18.27 


16.2 


2.5 


79±4 


2.2 


-0.07 


1261 


-17.38 


21.1 


2.5 


54±1 


2.5 


-0.27 


1308 


-15.41 


9.1 


2.0 


47±1 


3.6 


-0.46 


1407" 


-15.64 


8.9 


8.9 


36±5 


5.0 


-0.51 


1422'' 


-17.07 


19.3 


1.5 


44±7 


7.1 


-0.48 


1491" 


-15.44 


10.2 


2.5 


47±4 


5.6 


-0.44 


1514 


-15.58 


13.7 


2.0 


42±3 


4.5 


-0.77 


1545'' 


-15.90 


11.2 


1.5 


59±5 


4.3 


-0.05 


1743 


-15.20 


10.4 


10.4 


44±8 


14.0 


-1.01 


1857 


-15.85 


20.8 


20.8 


46±16 


2.7 


-0.66 


1871'' 


-16.54 


6.9 


1.5 


68±3 


4.9 


-1-0.03 


2019 


-16.25 


15.4 


1.5 


44±1 


2.3 


-0.28 


2048" 


-16.85 


12.6 


2.5 


33±2 


3.9 


-0.44 


2050 


-15.55 


10.8 


2.5 


32±2 


2.6 


-0.54 



" Data from the HyperLEDA FITS archive (OHP CARELEC). 
'' IFU spectroscopic data (MPFS). 

and the Hyper Leda FITS Archive, the agreement between 
the derived kinematical and stellar population profiles is 
fairly good (see details in the following subsection on in- 
dividual galaxies). 

We quantitatively compar e our kinematical p rofiles 
with the results published in ISimien fc Prugn"ieil (|2002l ) 
which are available in the electronic form. In all the cases 
both, the measurements and the derived uncertainties are 
absolutely consistent between the two studies although they 
were obtained using completely different data analysis tech- 
niques. We notice that in order to derive the stellar popu- 
lation properties we need higher signal-to-noise ratio than 
required only to constraint kinematics, therefore we had 
to apply coarser binning to the data degrading the spa- 
tial resolution to a high er degree than that presented in 
ISimien fc Prugniell (|2002l ). 

Since no kinematical profiles are provided in the 



0.08 


0.03 


0.14±0.05 


4.8 


8.5±2.0 


-0.81±0.04 


0.09 


0.07 


...±... 


2.8 


3.6±0.4 


-0.27±0.03 


0.10 


0.04 


...±... 


3.6 


3.6±0.4 


-0.27±0.03 


0.04 


0.02 


0.06±0.05 


3.2 






0.23 


0.09 


...±... 


3.9 


2.4±0.3 


-0.15±0.03 


0.13 


0.04 


...±... 


2.4 


2.4±0.3 


-0.15±0.03 


0.06 


0.01 


0.09±0.06 


3.2 


4.9±0.7 


-0.30±0.03 


0.14 


0.04 


...±... 


3.6 


4.9±0.7 


-0.30±0.03 


0.27 


0.07 


...±... 


3.2 


9.3±2.9 


-0.79±0.05 


0.13 


0.07 


...±... 


3.1 


2.5±0.3 


-0.17±0.03 


0.10 
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computer-readable form in the other publications presenting 
dE/dSO kinematics, we were able to make only qualitative 
comparison with the published kinematical profiles for some 
objects. 

Ivan Zee et all (|2004bh present two kinematical profiles 
for every galaxy observed in the course of the Palomar 
dE project, corresponding to the spectra obtained in the 
blue and red arms of the Double Spectrograph. The red 
arm data, having higher spectral resolution and sharper and 
deeper absorption lines of the Ca II triplet result in higher 
quality kinem a tical d ata, which is clearly seen in Fig 3 of 
Ivan Zee et al.l (|2004bt ). It is also clear from the same fig- 
ure that the velocity dispersion values obtained from the 
blue arm spectra are generally systematically above those 
derived from the red arm data, which is probably caused 
by the imperfections of the data analysis in low-cr regime, 
where measured velocity dispersions are similar to or lower 
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Table 3. Kinematical and stellar population properties of dwarf early-type galaxies: external regions. The columns are as follows: object 
name, velocity dispersion in the periphery, age, metallicity, 2 columns representing uncertainties of stellar population parameters (see 
text), [Mg/Fc] abundance ratio, B-band mass-to-light ratio, metallicity gradients per re and per kpc where positive values denote the 
decrease of [Fe/H] outwards. 
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" Data from the HyperLEDA FITS archive (OHP CARELEC). 
IFU spectroscopic data (MPFS). 

than the instrumental resolution. Our analysis of the blue 
arm spectra in terms of kinemati cs turns to be much m ore 
precise than what is presented in Ivan Zee et al.1 (|2004bt ). In 
the central regions of the galaxies our velocity dispersion 
profiles agree remarkably well with those derived f rom th e 
higher resolution red arm spectra by Ivan Zee et al.1 (|2004bD . 
whereas in the peripheral parts we usually reach even higher 
quality of measurements evident from the lower scatter of 
the values. The two exceptions are VCC 956 and VCC 2050, 
where velocity dispersions are too low to be extracted from 
blue arm spectra. 

The kinematical profiles of several galaxies included in 
our sample were presented in the literature. The rotation 
profiles ag ree weU for VCC 54 3 and VCC 1010 displayed 
in Fig 2 of IPedraz eraP lj2002h . while in VCC 1122 we see 
somewhat higher degree of rotation. The velocity dispersion 
values are consistent for all three objects. Their data were 



obtained with the IDS spectrograph at the 2.5 Isaac Newton 
Telescope and had ~2.5A FWHM spectral resolution in the 
blue spectral range, which is very similar to the datasets we 
discuss in this paper. 

Seven galaxies included in our sample, VCC 543, 856, 
917, 1036, 1087, 1261, and 1308 were observed with the 
Keck ESI sp ectrograph, and the kin ematical analysis was 
presented in ICeha et al.1 (|2002l . bOOSj ). The data had con- 
siderably higher spectral resolution than ours {R ~ 10000) 
but were localised in the central regions of the galaxies (out 
to 10 arcsec from the centres) due to a relatively short slit. 
Both, radial velocity and velocity dispersion profiles for the 
seven galaxies generally agree, however we notice that the 
specific details (e.g. cr-drops) look smeared in the ESI data 
compared to the Palomar DS spectra, although one would 
expect the opposite, given better atmosphere conditions at 
Mauna Kea. Presently, we do not have an explanation to 
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Figure 3. Comparison of age (top panel) and metallieity (bottom 
panel) measurements in the central parts of the galaxies made 
on the SDSS DR6 spectra and on the long-slit data. Different 
symbols correspond to the data from the Palomar dE project 
and HyperLeda FITS Archive. 



this effect, but our guess is that it may be connected to 
the atmosphere dispersion which would be an order of a 
few arcsec on a very wide wavelength range of ESI even at 
moderate air-masses and might smooth the kinematical de- 
tails if one fits the spectra in the entire wavelength range 
of the instrument. There is no mentioning about the at- 
mosphere dispersion corrector in the ESI manual available 
online, nei ther any info rmation about how it was taken into 
account in iGeha et all (t2002, . .2003 ). Comparing the shapes 
of the velocity di spersion profiles fo r VCC 917, 1036, and 
12 61 in Fig 5 in ICeha et aD lj2002h and Fig 2 and Fig 3 
m iGeha et al.l (j2003h . to the data presented in our study 
one would understand why Geha et al. did not mention the 



kinematically-decoupled cores in these galaxies: they were 
smoothed out. 

No spatially-resolved information about the stellar pop- 
ulation of the galaxies in our sample is available in the 
literature. Qualitative comparison of th e int egrated values 
with th ose presented in lCeha et al.l (|2003l ) and lvan Zee et al] 
(|2004al ) based on the analysis of the line-strength indices re- 
sults in a good agreement. The statistical comparison of the 
dE/dSO stellar population properties with another samples 
will be given in one of the next papers of the series. 



4.2 Notes on Individual Galaxies 

Here we briefly describe the results obtained for every 
galaxy. All cases with spatially-resolved information avail- 
able are illustrated with the profiles or two-dimensional 
maps of kinematical and stellar population properties, ex- 
cept VCC 1871, where th e colour figures are available in 
IChihngarian et all (j2007ch . 

VCC 178 = IC 3081 (Fig©. This flattened nucleated dwarf 
galaxy exhibits significant rotation and quite flat velocity 
dispersion profile inside 1 re- No em bedded substructures 
were revealed by iLisker et al.1 (|2006l ): VCC 178 is listed in 
Appendix E there. There is a significant metallieity gradient. 
Age distribution also demonstrates smooth changes along 
the radius, being younger in the centre (7 Gyr) than in the 
periphery (11.5 Gyr). 

VCC 389 = IC 781 (Fig[5|. This n ucleated dE palaxy host- 
ing a ■ 



"probable disc" reported by iLisker et al 



observed twice with the CARELEC spectrograph at slightly 
different positional angles. Both datasets have quite low 
signal-to-noise ratios, however the derived profiles of the 
stellar population parameters along the radius agree well. 
We detected neither statistically significant difference of the 
stellar population properties in the nucleus and in the main 
galactic body, nor metallieity gradient. 
VCC 437 = UCC 7399A (Fig[6l). A very steep metallieity 
gradient is observed beyond 4 arcsec from the centre. Metal- 
lieity at 1 re is almost 1 dex lower than at 4 arcsec. Inter- 
estingly, at ^3 arcsec the q' — z' c olour profile of the galaxy 
presented in lFerrarese et all (|2006l ) changes its gradient, and 
the ellipticity changes drastically. This suggests the presence 
of an embedded structure in the central part, however, not 
clearly detected on the unsharp-masked SDSS images, which 
can be explained if it has nearly face-on orientation. This ob- 
ject is one of the cases, where the luminosity-weighted age 
also changes along the radius: the population becomes older 
at larger radii, creating a "conspiracy" effect in the broad- 
band optical colours, making colour profiles looking nearly 
flat. 

VCC 490 = IC 783 (Fig [T]). This nucleated dwarf 
galaxy exhibits prominent sp iral arms in the stellar disc 
(jBarazza et"all |2002| : ILisker ot al.. 200^ ). The HFA data 
available for this object have very low signal-to-noise ratio 
outside the nuclear region, therefore the stellar population 
properties cannot be derived in the periphery of the galaxy. 
Deep IFU observations were carried out with the MPFS 
spec trograph in 2004 and prese nted in IChilingarianI (|2006l ) 
and IChilingarian et all (|2007cl ). The luminosity- weighted 
age in the nuclear region of VCC 490 is as young as 3 Gyr, 
whereas the main galactic body is older than 10 Gyr. There 
is a hint of a large-scale solid body rotation of the spiral disc 
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Figure 4. Kinematics and stellar populations of VCC 178 (Palo- 
mar DS). Four panels (from top to bottom) present: radial ve- 
locity, velocity dispersion, SSP-equivalent age, and metallicity. 
Results from the fitting of a SDSS spectrum are shown by filled 
gray circle. Positional angle of the slit and target signal-to-noise 
ratio used to adaptivcly bin the data are indicated. Dotted verti- 
cal linos indicate the photometric centre and ±re. 



llSimien fc Prugnie]|[2003) . The galaxy resides in the periph- 
ery of the Virgo cluster, close in projection to the luminous 
spiral M 100. Their radial velocities are different by only 
270 km s"^ which may suggest the belonging of VCC 490 
to the M 100 group. 

VCC 543 = UGC 7436 (FigEJ. This flattened rotationally 
supported galaxy was a subject to the several kinematical 
studies referenced above. It is one of the three galaxies with 
the data available both in the Palomar dE project and in 
the HFA. The kinematical and stellar population profiles are 
consistent between the datasets. The age profile in this ob- 
ject is nearly fiat with a mean value of ~5.5 Gyr, while the 
metallicity exhibits a slight gradient . No embedded str uc- 
tures were detected in this object bv lLisker et~aLl (|2006h . 
yCC 545 = IC 1 83 A No substructures were detected by 
iLisker et al.l (|2006l ) in this faint nucleated dE galaxy. The 
signal-to-noise ratio of the HFA data was very low so we were 
able to determine only integrated properties of the stellar 
populations and mean velocity dispersion, which agree well 
with the values derived from fitting the SDSS DR6 spec- 
trum. The galaxy is located close in projection to VCC 490 
and M 100 and its radial velocity is quite close to that of 
VCC 490. 

VCC 634 = NGC 4328 (Fig©. This relatively bright ro- 
tationally supported nuclea ted dwarf galaxy does not ex- 
hibit embedded structures (jLisker et al.l 1200a ). We see the 



Figure 5. Kinematics and stellar populations of VCC 389 (OHP 
CARELEC). See Fig |4] for details. Data from two datasets for 
slightly different positional angles are shown with diamonds and 
triangles. 



intermediate-age (3 Gyr) quite metal-rich (—0.2 dex) stel- 
lar population without significant changes of the properties 
within 1 re. Low signal-to-noise ratio of the data did not al- 
low us to measure the metallicity behaviour at larger radii. 
VCC 856 = IC 3328 (Fig \W^. This was the first dwarf 
dE/dSO g alaxy, where low-c ontrast spiral arms were de- 
tected by IJerien et al.] (|2000l) . The data allowed us to de- 
rive the stellar population properties out to 1 r^. The age 
profile is fiat with the mean value of about 4 Gyr, whereas 
the metallicity distribution exhibits a gradient. With the 
presently available data we were not able to detect any sta- 
tistically significant difference between the ages of the stellar 
populations in the nuclear region and in the spiral disc. 
VCC 917 ^ IC 3344 ( FigHTj) . No substructures were re- 
vealed by iLisker et all (|2006l ') in VCC 917. This fiattened 
object with a very little rotation was considered as one of 
the examples of the galaxies supported by the anisotropic ve- 
locity dispersions. We see a kinematically decoupled central 
component with the radius of about 3 arcsec (~0.25 kpc) 
rotating in the other sense compared to the main galactic 
body. It is associated with the drop in the velocity disper- 
sion profile from ~40 to ~20 km s^^ and younger and much 
more metal-rich stellar population compared to the periph- 
ery. This is the first example of the KDC in a dwarf galaxy 
supported by anisotropic velocity dispersions. 
VCC 965 = IC 3363 (Fig [T ^. No substructures are re- 
ported bv lLisker et al. 1 l|2006l ) in this faint fiattened nucle- 
ated galaxy. This is one of a few cases, where the internal 
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Figure 6. Kinematics and stellar populations of VCC 437 (Palo- 
mar DS). See Fig |4] for details. 



velocity dispersion is too low to be measured. The stellar 
population in the nuclear region turns to be very diflerent 
from the periphery of the galaxy both in terms of young age 
and high metallicity. The velocity dispersion profile of this 
o bject derived from the red arm DS spectra and presented 
in lvan Zee et al.1 (|2004bh demonstrates a central velocity dis- 
persion drop. 

VCC 990 = IC 3369 (Fig[I31). iLisker et al.l (|2006D report 
the detection of an inclined disc in this flattened nucleated 
galaxy. Inside 1 re we see a quite strong metallicity gradi- 
ent from about —0.3 dex to —0.8 dex while the age remains 
nearly constant at a level of ~5 Gyr. The velocity disper- 
sion distribution exhibits a local minimum in the centre as- 
sociated with a peculiarity on the radial velocity profile. We 
detected faint emission lines (H/3, [OIII]) in the fitting resid- 
uals suggesting the presence of ionised gas in this galaxy. 
VCC 1010 = NCC 4431 (Fig[l4l). Bar and disc were de- 
tected in this flat tened rotationally supported galaxy by 
iLisker et al.l (|2006l ). The data available for this object in the 
HFA allowed us to trace the behaviour of the stellar pop- 
ulation parameters out to re/2. Both, age and metallicity 
distribution look flat showing no signiflcant changes along 
the radius. 

VCC 1036 = NCC 4436 (FigfTS)) . The inclined disc is re- 
ported by ILisker et al.l (|2006l ) . Nearly solid-body rotation 
out to 1 re changes into the flattening of the rotation curve 
outwards. The velocity dispersion stays constant beyond 
0.3 Te at a level of ~50 km s~^ . This galaxy exhibits a promi- 
nent KDC. It is one of the three objects i ncluded both in the 
Palomar dE sample and in the sample of ISimien fc Prugniell 






Figure 7. Kinematics and stellar populations of VCC 490 
(MPFS). First row: radial velocity and velocity dispersion, sec- 
ond row: age and metallicity. Adaptive binning with the target 
signal-to-noise ratio of 15 was applied. 
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Figure 8. Kinematics and stellar populations of VCC 543 (Palo- 
mar DS and OHP CARELEC). Two curves on each panel cor- 
respond to the two different datasets. The target signal-to-noise 
ratios were 20 and 10 for the Palomar DS and OHP datasets 
respectively. See Fig |4] for details. 
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Figure 9. Kinematics and stellar populations of VCC 634 (OHP 
CARELEC). See Fig H] for details. 



Figure 11. Kinematics and stellar populations of VCC 917 (Palo- 
mar DS). Sec Fig |4] for details. 
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Figure 10. Kinematics and stellar populations of VCC 856 (OHP 
CARELEC). See Fig H] for details. 



Figure 12. Kinematics and stellar populations of VCC 965 (Palo- 
mar DS). See Fig |4] for details. 
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Figure 13. Kinematics and stellar populations of VCC 990 (Palo- 
mar DS). See Fig |1] for details. 
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Figure 14. Kinematics and stellar populations of VCC 1010 
(OHP CARELEC). See FigHfor details. 



Figure 15. Kinematics and stellar populations of VCC 1036 
(Palomar DS and OHP CARELEC). The target signal-to-noise 
ratios were 30 and 10 for the Palomar DS and OHP datasets 
respectively. See Fig [S] for details. 



(|2002l ). Both datasets agree well on the presence of the kine- 
matically decoupled central component. However, due to 
worse seeing conditions it is not seen as counter rotation 
in the HFA data, but as a flat region of the radial velocity 
proflle. A slight depression is seen in the velocity dispersion 
proflle in the region corresponding to the KDC. It corre- 
sponds to a pronounced bump in the metallicity profile with 
the values reaching almost [Fe/H]=0.0 in the SDSS data. 
SSP-equivalent age (~3Gyr) in this part is ~2 Gyr younger 
than in the surrounding regions. The metallicity keeps con- 
stant at a level about —0.3 dex out to 8 arcsec, then it starts 
to decrease steeply. This is the only object, where the cen- 
tral value of the a/Fe abundance ratio (—0.05 dex) differs 
strongly from the values in the periphery (-1-0.17 dex). 
VCC 1075 = IC 3383 (Fig[l6l). This faint flattened nucle- 
ated dwarf galaxy harb o urs n o embedded substructures ac- 
cor ding to iLisker etHI (120061 ). The signal-to-noise ratio of 
the data is quite low, making difficult to assess the variations 
of the stellar population parameters along the radius. How- 
ever, at the present level of precision we do not detect any 
statistically significant gradients of age neither metallicity. 
VCC 1087 = IC 3381 (Fig [HI). This flattened nucleated 
dwarf with no substructures reported bv lLisker" et all (j2006[ ) 
has very little if any rotation. The low signal-to-noise ratio 
of the data makes difficult to study the variations of the stel- 
lar population along the radius. We notice that the stellar 
population parameters derived from the SDSS DR6 spec- 
trum fitting are very different from those estimated from 
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Figure 16. Kinematics and stellar populations of VCC 1075 
(Palomar DS). See Fig |4] for details. 



Figure 17. Kinematics and stellar populations of VCC 1087 
(OHP CARELEC). See Fig H for details. 



the HFA spectra. This, however, may be explained as the 
seeing effect, which smoothed the nucleus. 
VCC 1122 = IC 3393 (Fig [TH]). This strongly flattened 
galaxy with the insufficient degree of rotation to be consid- 
ered as rotationally supp orted does not exhi bit embedded 
substructures according to lLisker et al.l (|2006h . We have the 
data from both Palomar dE project and the HFA. The kine- 
matical profiles agree quite well. There is a peculiarity seen 
in the radial velocity profile 2 arcsec south-east of the nu- 
cleus, however presently we do not have enough evidences 
to interpret it as a kinematically decoupled component. The 
stellar population properties were determined out to 2 
from the Palomar dE data. The nuclear region of VCC 1122 
is clearly distinct in both age and metallicity profiles. We 
also see a remarkable metallicity gradient, which is, however 
still insufficient to explain the high metallicity in the nuclear 
region. The age stays constant along the radius beyond re/2. 
There are faint H/3 and [0111] emission lines revealed in the 
fitting residuals localised in the central part of the galaxy. 
VCC 1183 = IC 3413 (Fig[T9l). Thi s fiattened n u cleate d dE 
galaxy harbours a bar reported bv lLisker et al.l (|2006l ). We 
observe some rotation. Low signal-to-noise ratio of the data 
did not allow the detailed studies of the stellar populations, 
however we notice the intermediate age (~4 Gyr) and rela- 
tively high metallicity (—0.3 dex) without evident changes 
along the radius. 

VCC 1250 = NCC 4476 (Fig This galaxy hosts a 
large embedded disc with the spiral arms, dust lanes, and 
HII regio ns clearly visibl e on t he HST ACS images pre- 
sented in iFerrarese et"al] (120061 ). Our data covers the re- 



gion completely dominated by the disc. We see prominent 
emission lines suggesting on-going star formation. The lumi- 
nosity weighted age is about 2 Gyr and the metallicity is as 
high as solar without statistically significant changes along 
the radius. The _B-band absolute magnitude of VCC 1250, 
— 18.27 mag places it at the limit of the informal dwarf/giant 
separation. However, it still can be classified as a low- 
luminosity early- type galaxy, that's why we included it into 
our sample. 

VCC 1261 = NCC 4482 {Fig \2T]i . No substructures were 
reported bv lLisker et al.l (|2006l ) in this galaxy. This object 
looks like an enlarged version of VCC 917, giving another 
example of a small rapidly rotating KDC in a fiattened 
spheroid with a very low, if any, large-scale rotation. As in 
VCC 917, the KDC is associated with a velocity dispersion 
drop, younger and much more metal rich stellar population 
than in the peripheral parts of the galaxy. Faint emission 
lines are seen in the fitting residuals, suggesting the pres- 
ence of warm ISM in this object. 

VCC 1308 = IC 3 437 (Fig[22l) This galaxy was classified as 
non- rotating " bvlGeha et all (j2003lV but as "rotating" by 
Ivan Zee et all (l2q04bl). No substructures were reported in 
it bv lLisker et al.l (|20od ). We indeed see some rotation, al- 
though the radial velocity profile is quite peculiar in the cen- 
tre, as well as a remarkable metallicity gradient. The stellar 
population in the central region of the galaxy looks younger 
then in the outer regions along the radius. 

VCC 1407 = IC 346I Very low signal-to-noise ratio of the 
data does not allow to study the variations of the stellar 
population properties along the radius. However, we stress 
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Figure 18. Kinematics and stellar populations of VCC 1122 
(Palomar DS and OHP CARELEC). The target signal-to-noise 
ratios were 30 and 10 for the Palomar DS and OHP datasets 
respectively. See Fig [S] for details. 



a very good agreement of the age and metallicity estimates 
for this object obtained from the fitting of the HFA and 
SDSS DR6 spect ra. No embedded s tructures were reported 
in this galaxy by lLiske~ et all (|2006l ). 

VCC 1422 = IC 3468 (Fig [231). An embedd ed structure is 
detect ed in this galaxy interpr eted as a bar bylBarazza et al.l 
|2002) or as an edge-on disc by lLisker et"al] (|2006l '). The IFU 
spectroscopy reveals complex kinematics with evidences for 
rotation only beyond 5 arcsec from the centre. There is a 
hint that stellar population in the region of the disc/bar is 
somewhat younger than in the periphery of the galaxy (a = 
54± 4 km s -\ t = 7.2 Gyr, [Fe/H]=-0.54 dex), although 
the difference is on the limit of detection. VCC 1422 is the 
only galaxy in our sample, where the very central region, 
dominated by the light of its nucleus exhibits older and more 
metal-poor stellar population than one of the surrounding 
disc/bar clearly visible on the maps. Its age and metallicity 
are found to be very close to those in the periphery of the 
galaxy. 

VCC 149 1 = IC 3486 (Fig|24l) . No substructures were re- 
vealed by iLisker et all (12009 )" in this non-nucleated galaxy 
demonstrating quite high degree of rotation for its moder- 
ate flattening. The stellar population is old (~8 Gyr) and 
metal-poor (—0.6 dex) without noticeable gradients. 
VCC 1514 = PCC 41726 (Fig[25l) . Possible inclined disc 
was reported bv lLisker et al.l (|2006l ) in this strongly flattened 



Figure 19. Kinematics and stellar populations of VCC 1183 
(OHP CARELEC). See Fig H for details. 



VCC1250, P.A.=203 deg, S/N= 10 



E 60 
°^ 40 



: 

! 1 — ' 


> : 














; 


> 

. <> 






: 









-30 -20 



10 20 



position, arcsec 



gala^xy. We observe some rotation. The age and metallicity 
estimates do not exhibit statistically significant variations 



Figure 20. Kinematics and stellar populations of VCC 1250 
(OHP CARELEC). See Fig H for details. 
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Figure 21. Kinematics and stellar populations of VCC 1261 
(Palomar DS). See Fig |4] for details. 



VCC1308, P.A.= 66 deg, S/N= 20 
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Figure 22. Kinematics and stellar populations of VCC 1308 
(Palomar DS). See Fig |4] for details. 
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Figure 23. Kinematics and stellar populations of VCC 1422 
(MPFS). Panels are the same as in Fig [7] Adaptive binning with 
the target signal-to-noise ratio of 15 was applied. 

along the radius, although the signal-to-noise ratio of the 
data was rather low. 

VCC 1545 = IC 3509 (Fig[26|. This galaxy was among the 
first, where the young stellar population in the central region 
was reported bv lChilingarian et al.l ()2007dl ) . It demonstrates 
very complex kinematics: we see a minor-axis rotation (i.e. 
kinematical decoupling) in the circumnuclear region of the 
galaxy with the depression in the velocity dispersion distri- 
bution there, while there is a large-scale major-axis rota- 
tion as well. The stellar population in the nuclear region is 
younger and considerably more metal-rich than in the main 
galactic body. 

VCC 1743 = IC 3602 This is a faint flattened dwarf galaxy 
with the faint extended nucleus seen on the HST ACS im- 
ages but not on the g round-based d ata. T he galaxy is listed 
in Appendix E of of iLisker et al.l l|2006h as an object "in 
which substructure other than a disc was found" . The data 
available to us have very low signal-to-noise ratio due to low 
surface brightness of the galaxy, therefore we were not able 
to derive detailed profiles of the stellar population proper- 
ties. However, outside the nuclear region the age looks to be 
old (~10 Gyr) and the metallicity quite low (—0.6 dex). 

VCC 1857 = IC 3647 The faintest object in the Palomar 
dE sample does n ot exhibit and substructures reported by 
iLisker et all (|2006l ). The signal-to-noise is not sufficient even 
for the precise integrated measurements of the stellar popu- 
lations, however at 1-a confidence level the galaxy has and 
intermediate age not older than 7 Gyr. 

VCC 1871 = IC 3653 An inclined stellar disc was d iscov- 
ered in the velocity field bv lChilingarian et al.l (|2007cl ) with 
a counterpart in the colour map. The galaxy has somewhat 
high surface brightness, the age distribution is nearly uni- 
form, whereas there is a pronounced peak in the metallicity 
map with the the values slightly exceeding solar. 

VCC 2019 = IC 3735 (Fig [HI). A "possible inclined disc, 
maybe warped or distorted" is reported in this galaxy by 
ILisker et alj (|2006l ). The isophotes demonstrate the change 
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Figure 24. Kinematics and stellar populations of VCC 1491 
(OHP CARELEC). See Fig g] for details. 
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Figure 25. Kinematics and stellar populations of VCC 1514 
(Palomar DS). See Fig |4] for details. 



of ellipticity and positional ang le (jPerrarese et al.ll2006t ). We 
see a central peak in the metallicity distribution and the cor- 
responding drop in the velocity dispersion profile, while the 
quality of the age measurements is not sufficient to clearly 
say whether the nucleus harbours younger stellar population 
than the outer parts of the galaxy. A noticeable metallicity 
gradient is observed. We detected faint emission lines (H/3 
and [OIII]) in the fitting residuals. 

VCC 2048 = IC 3773 (Fig There is an evidence for 
a bar or edge-on embedded disc reported by iLisker et al.l 
(|2006l ). This nucleated galaxy is strongly flattened. Its inner 
region, corresponding to the disc demonstrates fast solid- 
body rotation and looks kinematically decoupled from the 
outer parts of the galaxy rota ting in the same sense, similar 
to the structures observed bv lDe Riicke et ahl |200J) in the 
two dE galaxies in groups. There is a pronounced cr-drop in 
the central region of the galaxy. Stellar population properties 
do not show statistically significant changes along the radius 
apart from the little region 4 arcsec North of the nucleus, 
where the metallicity rises by ~0.3 dex compared to the 
nearby regions. Presently we cannot explain this feature, it 
may be connected to the globular cluster which might have 
fallen into the slit. 

VCC 2050 = IC 3779 (Fig f29l). Probab le incli ned disc is 
reported in this galaxy by iLisker et al.l (|2006l ). The ob- 
served velocity dispersion in this object is too low to be 
precisely measured. There is a significant metallicity gradi- 
ent, whereas the signal is not sufficient to detect any changes 
in the age distribution. The luminosity- weighted age in the 
central part is quite young (~2.5 Gyr) and we see faint emis- 
sion lines in the fitting residuals. 
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Figure 27. Kinematics and stellar populations of VCC 2019 Figure 29. Kinematics and stellar populations of VCC 2050 

(Palomar DS). See Fig|4]for details. (Palomar DS). Sec Fig|4]for details. 
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Figure 28. Kinematics and stellar populations of VCC 2048 
(OHP CARELEC). See FigHfor details. 



5 SUMMARY AND CONCLUSIONS 

We have presented the first large dataset of spectroscopi- 
cally derived radial profiles of stellar population parameters 
for dwarf early-type galaxies made by re-analysing some of 
the published spectroscopic data. Thanks to the usage of the 
NBuRSTS full spectral fitting techni que we have improved 
the stellar kinematics compared to Ivan Zee et al.l (|2004bt l 
reaching comparable quality of m easurements to those pre- 
sented in lGeha et al] (|2002l . 120031 ') but generally going out to 
1 re or further from the centre. We have performed the stel- 
lar population analy sis of the spectra where o nly kinematical 
data were available jSimien fc Prugniel|[2003 ). Our measure- 
ments of the nuclear stellar populations are in a good agree- 
ment with the results obtained from the fitting of SDSS DR6 
spectra. 

Significant metallicity gradients are often observed in 
dE/dSO galaxies, although in some of the objects the metal- 
licity distribution is completely constant along the radius. 
There is a tendency to the dichotomy in the distribution 
of the metallicity gradients - metallicity profiles are either 
nearly fiat or steeply decreasing along the radius, reaching 
—0.9 dex in VCC 437. In none of the cases we see a sta- 
tistically significant positive gradient. Age distribution along 
the radius usually exhibits fiat behaviour. If any gradients in 
age exist, they cannot be clearly detected using the existing 
datasets. 

Levels of «-enhancement as derived from the Lick in- 
dices exhibit no statistically significant changes along the 
radius and between the nuclei and main discs/spheroids in 
all objects, but VCC 1036, where the value in the nuclear 
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region is nearly solar ([Mg/Fe]=— 0.05 dex), whereas it stays 
super-solar (+0.17 dex) along the radius out to 1.2 re- 

Stellar populations in the nuclear regions are usually 
more metal-rich than in the main galactic bodies, and the 
central values of the luminosity-weighted metallicity are 
higher than what one would expect from the simple in- 
terpolation of the metallicity gradients toward the cen- 
tres. A remarkable fraction of galaxies (11 or 27) demon- 
strate luminosity-weighted ages in the nuclear regions to be 
younger than those of the main galactic bodies, with the 
difference exceeding 7 Gyr in VCC 917 and VCC 490. 

These chemically and evolutionary decoupled nuclear 
regions are always associated with the drops in the velocity 
dispersion profiles. Only in those cases, where the velocity 
dispersions are too low to be measured with the present 
data, the a-drops are not observed. In these cases, the mass- 
to- light ratios change significantly along the radius, pointing 
out to unrighteous use of constant M /L dynamical models. 

In certain objects (VCC 990, 1122, 1250, 1261, 2019, 
and 2050) we detected faint emission lines (H/3, [OIII]) in 
the nuclear regions, associated with the chemically and evo- 
lutionary decoupled cores. 

We discovered four kinematically decoupled rotating 
cores in bright members of our sample (VCC 917, 1036, 
1261, and 1545), two of them (VCC 917 and 1261) being 
hosted by the galaxies with very low if any major axis rota- 
tion. In all cases KDCs are associated with young metal-rich 
stellar populations and velocity dispersion drops. 

Finally, we see remarkable resemblance between inter- 
nal structure of dE/dSOs and of intermediate-mass and gi- 
ant early-type galaxies. Some low-luminosity objects after 
detailed analysis look like "scaled down" versions of more 
luminous counterparts (e.g. VCC 917/1261 and NGC 5813, 
VCC 1545 and NGC 4365), suggesting possible similarities 
between their evolutionary paths. Some physical processes, 
e.g. major mergers and secular evolution, are usually consid- 
ered to rule the life of giant early-type galaxies. We have to 
revise their applicability to explain the dE/dSO evolution. 
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